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ABSTRACT 

It  is  ob8erve<1  that  Vorobyev's  methor^  of  moments,  because  it  is  a 
specialized  form  of  Galerkin's  method,  is  easily  accelerated  by  the  use  of 
Sloan's  trick  of  taking  first  iterates.  The  improvement,  which  is  nearly 
costless  since  it  requires  no  new  quantities,  is  illustrated  in  a  numerical 
example . 
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SIGNIFICANCE  AND  EXPLANATION 


Operator  equations  have  countless  realizations  in  miscellaneous 
engineering  and  scientific  applications.  The  language  of  functional  analysis 
provides  a  powerful  and  unifying  tool  for  the  study  of  phenomena  associated 
with  the  computational  solution  of  such  equations.  The  functioral  analytic 
approach  avoids  tedious  details  and  allows  the  analyst  to  focus  on  major 
issues.  This  methodology,  though  it  may  not  produce  actual  software,  can 
provide  invaluable  conceptual  understanding. 

We  deal  with  the  approximate  solution  of  the  equation 

(J  -  X)u  ■  V  , 

where  the  operators  J  -  K  and  J  are  assumed  invertible,  by  the  use  of 
Vorobyev's  method.  This  technique  uses  the  n  +  1  "moments" 

Zq  ■  j“^v,  Z£  »  1  <  i  <  n  , 

in  order  to  get  an  approximate  solution  of  the  form 

^n  “  ^0*0  +  +  •••  +  Cn-I^n-I  • 

Under  proper  conditions,  the  sequence  (y^)  converges  to  the  true  solution 
u*  -  (J  -  K)"V. 

Sloan  has  shown  that,  given  an  approximate  solution  u,,  to  the  equation 
the  first  Iterate 

^  -  J'’ku„  +  zq 

is  usually  more  accurate  than  Ujj  itself.  Specifically,  under  adequate 
conditions,  the  sequence  of  first  iterates  {u^}  converges  to  u*  faster 
than  the  sequence  {u^j} .  The  aim  of  our  note  is  to  report  the  simple  but 
useful  observation  that  Sloan's  iterative  refinement  is  applicable  to 
Vorobyev's  method  and  that  the  computation  of  the  first  iterate, 

^n  “  *0  ^0*1  *  *•*  ^n-1*n  ' 

involves  only  known  quantities.  A  numerical  example  involving  an  integral 
equation  Illustrates  this  nearly  costless  improvement. 


ITERATIVE  RBFIMEMEST  OF  THE  METHOD  OF  MOMENTS 
Gtorq*  Ml«l* 

Ii«t  E  and  F  danota  Hllbart  apacaa«  Glvan  boundad  linear  operators  J,K  i  E  >  F 
and  an  element  v  e  Ft  consider  the  split  equation  of  the  first  kind 

(J  -  K)u  -  V  .  (1) 

Aseuna  that  J  has  a  boundad  Inverse  on  F  and  that  E  Is  a  compact  operator.  Equation 
(1)  Is  equivalent  to  the  aquation  of  the  second  kind. 

u  a  Mu  V  t  (2 ) 

whore  M  •  J~^K  la  a  compact  operator  and  w  ~  J~^v.  The  aim  of  this  note  Is  to  report 
the  observation  that  Vorobyev's  method  of  moments  (8]  can  be  advantaqeously  Improved  by 
Sloan's  Iterative  refinement  [6].  Since  the  quantities  needed  for  the  latter  procedure 
are  conveniently  computed  In  the  former  method,  the  coat  of  the  Improvement  Is 
Inslqnlf leant.  After  a  brief  description  of  the  pertinent  procedures,  we  Illustrate  the 
effect  of  the  Iterative  refinement  on  a  Fredholm  Integral  equation  with  a  Green's  kernel. 

Sloan's  Iterative  Refinement 

Me  first  need  to  describe  Galerkln's  method  for  equation  (1).  Let  F^  •  *^0'  vhere 
Ejj  Is  a  flnlte-dlmenslonal  subspace  of  E,  and  let  be  the  orthogonal  projection  of 

F  onto  F„.  The  Galerkln  method  consists  of  solving  the  approximate  equation 

-  *’n’'>'*n  “  «n  '  ^n  • 

This  equation  Is  equivalent  to  the  equation 

"n  “  On”»n  ^n''  ' 

where  Pj^  «  la  the  orthogonal  projection  of  E  onto  Ej,.  A  standard  result 

*Addresss  P.  O.  Box  9208,  Marina  del  Rey,  CA  90295. 


stat«8  that  If  J  -  K  has  a  bounded  Inverse  and  0„x  *  x  for  every  x  e  E  as  n  ♦ 


then,  for  large  enough  n,  the  operator  I  -  Q„M  has  a  bounded  Inverse  on  E„t  the 
approximate  solutions  u„  -  (I  -  0„M)”'q„w  converge  to  the  true  solution 
u*  “  (J  -  K)”V,  and  the  error  bounds 


lu  -  u  I  <  I  (I  -  0  M)  'l  •  lu  -  Q  u 
n  n  Ti 


(4) 


are  valid.  We  refer  to  Krasnosel'skll  et  al.  (4]  for  a  proof. 

Sloan  [6]  has  shown  that  the  seguence  of  first  Iterates 

u^  «  Mu„  +  w  (5) 

converges  more  rapidly  to  the  true  solution.  Moreover,  once  Is  known,  the  iterate 

u^  can  be  computed  with  little  cost.  Indeed,  letting  e^,...,e,j  constitute  a  basis 

for  £^,  since  the  Galerkin  method  requires  the  quantities  In  the  setting-up 

n  n 

of  the  n  X  n  linear  system  needed  to  find  u^^  »  J  ^1*1'  iterate  5^  “  I  ^I'^i  "*■ 

1  1 

Is  easily  obtained. 

To  see  the  improvement  of  u^  over  u^,  we  proceed  as  follows.  Define  an  operator 
by 


«n  -  «Cn  • 

Then,  (3)  and  (S)  imply  that  Uj^  is  a  solution  of  the  equation 

%  -  «n^  *  ''  • 

Verify  now  that 

(I  -  M„)(u*  -  u^)  -  (M  -  M„)(u*  -  Q„u*)  . 

From  this  expression,  we  get  the  error  hound, 

»  _  _i  *  ♦ 

»u  -  u  I  <  I (I  -  M  5  >  •  *M  -  M  I  •  «u  -  0  u  I  , 

n  n  n  n 

which  shows  improvement  over  the  error  bound  in  (4).  See  Sloan  (6)  for  details. 

Recent  applications  of  the  iterative  refinement  to  Integral  equations  can  be  found  in 
Chandler  [1],  Craham  [3],  and  Spence  and  Thomas  [7).  A  functional  analytic  overview  of 
the  iterated  Galerkin  procedure  is  given  in  Chatelin  and  Lebbar  [2)  and  in  Schock  [5). 


/ 


Woroby«v»»  Method  of  MonentB 

He  turn  our  attontlon  back  to  aquation  (1)  antt  Its  aqulvalant  formulation  (2). 
Consider  the  n  1  ’moments” 

Sq  -  w,  <1  ~  ***  '  ^n  ~  **^-1  * 

The  usual  Neumann  Iterates  »  ^n-1  *  ttartinq  mlth  Xq  •  tt,  are  qlven  by 

*n  “  *0  +  *1  +  ♦  *n  • 

In  Vorobyev's  method,  one  uses  the  n  >  1  moments  (6)  In  order  to  find  a  linear 
combination  of  the  first  n  moments. 

^n  “  5o*0  +  5l*1  *  5n-1*n-1  • 

The  algorithm  proceeds  like  this: 

1.  Set-up  and  eolve  the  n  ”  n  linear  system 

n-1 

J  <*£,*. >o.  ■  -<*£,*„>.  0  <  1  <  n  -  1  ,  (8) 

1-0 

for  ^o''*''^n— 1' 

2.  evaluate  recursively  the  coefficients 

Cq  —  1  —  cOq.  —  o®^.  I^l^n  —  1. 

where  c  -  VCOq  +  ®n-1 

3.  evaluate  y^  using  (7). 

It  turns  out  that  the  algorithm  Is  equivalent  to  Galerkln’s  method 

Vn  "  On’*yn  *  '*•  ^n  *  ^n  '  *** 

where  is  the  orthogonal  projection  of  E  onto  E„  “  span(zg. .  •  •  rS„_^ )  •  'Hio  elements 

of  E.,  are  of  the  form  y„  -  q(M)zg.  where 

q(t)  •  Cg  +  C^t  +  •••  ♦  e„_it""’  (tO) 

is  an  arbitrary  polynomial  of  degree  <  n  -  1.  “nio  residual  la  given  by 

*0  ■  Vn  ’’yn  *  P<")*o  ' 

where 

p(t)  -  1  -  (1  -  t)q{t)  -  cCOg  ♦  o,t  ♦  •••  +  t")  .  (11) 


-y- 


-3- 


To  that  tha  approxlmata  oparator  aquation  (9)  raprasanta  Vorobyav'a  mathod,  taka 
t  "  1  In  (11)>  tharaby  showing  that  c  ■  ^/(*o  **•  ®n-1  ■*■')»  than  substltuta 

(10)  and  (11)  and  cos^re  coafflclanta  of  Ilka  powara  of  t  to  gat  tha  racurranca 


^  ralatlon  for 

j  If  tha  alamanta  ara  llnaarly  Indapandant  than  tha  llnaar  systan  (8)  la 

unlqualy  solvabla.  Otharwlsa>  tha  oparators  M  and  eolnelda  on  tha  subspaea 

and  thara  la  than  no  gain  to  ba  aada  by  using  tha  aathod.  If  tha  oparator  J  -  K  has  a 
boundad  invarse  than,  for  sufflclantly  larga  n,  aquation  (9)  has  a  uniqua  solution  y^ 

I  and  tha  saquanca  {y^}  convargas  to  tha  solution  (J  -  X)~V  fastar  than  any  gaonatric 

prograsslon.  Wa  rafar  to  Vorohyav  (8)  for  datails*  Mlscallanaous  applications  of  tha 
method  can  ba  found  In  Chaptars  V  and  VII  of  tha  citad  taxt. 

i 

Bafinamant  of  tha  Method  of  Moment 

Since  Vorobyev's  method  la  a  Galarkln  method  on  tha  subspaea 
E„  ••  span(Zg,MsQ, . . .  ,N"~^Zg) ,  Sloan's  Iterative  procedure  can  ba  applied.  In  view  of  (6) 

and  (7),  tha  desired  Iterate  Is  given  by 

Yn  *  5o*1  ^1*2  **•  5n-1*n  *  •  * 

Tha  refinement  Is  particularly  advantageous  since  tha  computation  of  (12)  requires  only 
already  known  quantities. 

Example .  Let  E  •  F  <•  L  [0,1]  and  J  -  I.  Consider  tha  Fredholm  Integral  aquation 
of  tha  second  kind 

u(a)  -  Ku(a)  +  e/2  , 

1 

K  J  E  ♦  E,  Ku(8)  -  /  k(s,t)u(t)dt  , 

0 


where  the  kernel  Is  given  by 


k(s,t) 


J  (2  -  t)e,  e  <  t 


B 


Til*  exact  solution  is  u  (s)  •  sin  apply  Vorobyev's  method  with  three  moments 

(n  -  2)t 


We  find  that: 


*0^*^  “  "2  ' 

z^(n)  -  Kzg(s)  “  ' 

8^5  3 

Z2(s)  -  K*,(s)  -  —  (3s=>  -  20s’  +  31s) 


<*0**0*  “  "12  * 


8  7 

<*0**1>  -  12  •  30 


8*  71  8^  517 

<*0**2>  -  <*1**1>  -  ^  *  mo  '  -  12  '  37800  ' 


a„  >2e.a  --LL82 

0  30  1  1260 


■  216  _  517-3 

7e»o  ^  -«r“i  -  mo  ® 


a8  -58 

‘O  ■  945  '  “l  "  18  ' 


r  7560  -  525V  .  _ 7560 _ 

■0“  2  4''l"  2  4 

7560  -  525ir<  +  4ir  7560  -  525*^  +  4*' 


Consequently,  the  desired  approxltaatlon  Is 


y2(e)  -  CoZo(*)  5i«i(*)  -  as  -  bs’  , 


Sloan's  refinement. 


•  ■'•552749  , 

b  -  -jg-  "  0.561564  . 


y2(s)  -  1Cy2(s)  +  v(s) 


requires  the  quantities  Kzg  •  and  Kz^  •  Z2  which  are  already  known.  He  find  that 

y2(e)  -  Cgz^(s)  tCi*2**'  +  v(s) 

«  cs  -  ds^  ♦  es5  , 


-•'V-  .• 


■5- 


c 


d 


e 


^^-n520 


48 

li’' 

3840 


576 


0.6385423 


r  •  0.06828021  . 


TABLE 


s 

u  (S) 

y2(«) 

y2<*> 

u  -  y2 

u  -  y2 

0.00 

0.00000 

0.00000 

0.00000 

0.00000 

0. 00000 

0.20 

0.30902 

0.30606 

0.30876 

0.00296 

0.00026 

0.40 

0.58779 

0.58516 

0.58754 

0.00263 

0.00025 

0.60 

0.80902 

0.81035 

0.80901 

-0.00133 

0.00001 

0.80 

0.95106 

0.95468 

0.95116 

-0.00362 

-0.00010 

1.00 

1.00000 

0.99118 

0.99998 

0.00882 

0.00002 

The  above  table  Illustrates  the  nearly  free  improvement  of  the  iterate  yj  over  y2 
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